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Abstract 

We predict and study theoretically a new nonlinear electromagnetic phenomenon in a sample of 
layered superconductor of finite length placed in a waveguide with ideal walls. Two geometries are 
considered here: when the superconducting layers are parallel or perpendicular to the waveguide 
axis. We show that the transmittance of the superconductor slab can vary over a wide range, from 
nearly zero to one, when changing the amplitude of the incident wave. Thus, one can induce the 
total transmission or reflection of the incident wave by just changing its amplitude. Moreover, the 
dependence of the superconductor transmittance on the incident wave amplitude has a hysteretic 
behavior with jumps. The considered phenomenon of self-induced transparency can be observed 
even at small amplitudes, if the wave frequency u is close to the cutoff frequency for linear Josephson 
plasma waves. 

PACS numbers: 74.78.Fk, 74.50. +r, 74.72.-h 
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I. INTRODUCTION 



There is considerable interest in studies of metamaterials and nanostructures with 
unusual electromagnetic properties. Layered high-temperature superconductors, e.g., 
B^S^CaC^Os+s single crystals, definitely belong to this type of materials. Experimental 
studies (see, e.g., Refs. Hp]) have shown that the electrodynamics of layered superconductors 
can be described by a theoretical model, which assumes that thin superconducting CuC>2 
layers (with a thickness s of about 2-3 A) are coupled through thicker dielectric layers (with 
a thickness d of about 15 A and a dielectric constant e ~ 15) via an intrinsic Josephson 
effect. Due to the layered structure of P^S^CaC^Og^ and similar compounds, they sup- 
port the propagation of specific electromagnetic waves, the so-called Josephson plasma waves 
(JPWs) (see, e.g., Refs. [3|H and references therein). These waves belong to the terahertz 
frequency range, which is very important for various applications but not easily accessible 
with modern electronic and optical devices. This technological perspective provides a strong 
motivation for research on these waves. From a scientific point of view, the interest in lay- 
ered superconductors is mainly related to the specific type of plasma formed, the so-called 
Josephson plasma. An essential property of the Josephson plasma is the strong anisotropy 
of its current-carrying capability. This anisotropy manifests not only in the difference of the 
absolute values of the flowing currents (currents along the crystallographic ab-plane are two 
orders of magnitude higher than those along the c-axis) but also in their physical nature. 
Indeed, the current along the layers has the same nature as in the usual bulk superconduc- 
tors and can be described within the framework of London's theory, while the current along 
the c-axis is of Josephson origin. 

A. Waves in usual plasmas versus Josephson plasma waves 

In the Josephson plasma, not only phenomena common to other types of plasmas can be 
observed but also those specific for layered superconductors. As in usual plasmas, there is a 
gap in the spectrum of Josephson waves. JPWs can only propagate with frequencies higher 
than the threshold Josephson plasma frequency uj. As was theoretically demonstrated in 
Refs. I5]lt3| the surface Josephson plasma waves (SJPWs) can propagate along the interface 
between layered superconductors and vacuum, similarly to usual plasmas. The excitation 
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of these waves leads to various resonant phenomena^ similar to the Wood anomalies well 
known in optics (see, e.g., Refs. I8 HTU1) . However, contrary to usual plasmas, SJPWs can 
propagate with frequencies not only below the plasma frequency but also above itP. The 
Josephson plasma can also exhibit properties characteristic for left-handed media: a negative 
refractive index for terahertz waves can be observed at its interface with vacuumP^. As 
was shown in Ref. [121 phenomena similar to Anderson localization and the formation of a 
transparency window for terahertz waves can be observed in layered superconductors with 
a randomly-fluctuating value of the maximum Josephson current. 

B. Nonlinearities 

Since the current along the c-axis is of Josephson origin, the electrodynamic equations for 
layered superconductors are nonlinear. This can lead to a number of non-trivial nonlinear 
effects accompanying the propagation of JPWs, e.g., slowing down of light 13 , self-focusing 
of terahertz pulses^^, excitation of nonlinear waveguide modes^, as well as self-induced 
transparency of the slabs of layered superconductors and hysteretic jumps in the dependence 
of the slab transparency on the wave amplitude^. 

C. Finite-size samples 

It should be noted that the majority of the theoretical studies of JPWs consider infinite 
samples. However, the sizes of layered superconducting samples used in experiments are 
comparable or even smaller than the wavelength of the terahertz radiation. Obviously, in 
this case the sample cannot be treated as infinitely large. This means that the theory should 
consider finite sample sizes. 

In this paper we discuss a nonlinear phenomenon that can be observed in a finite slab of 
layered superconductor placed inside of a rectangular waveguide with ideal walls. We have 
considered two geometries: when waves are propagating across the layers or along them. It 
turns out that the transmittance of the slab, being exposed by an incident wave from one 
of its sides, depends not only on the wave frequency, but also on the wave amplitude for 
both of the geometries considered here. Therefore, a slab of fixed length can be completely 
transparent (when neglecting dissipation) for waves of certain amplitudes, and nearly totally 
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reflecting for other amplitudes. Moreover, the dependence of the transmittance on the wave 
amplitude shows hysteretic behavior with jumps. The transmittance of the slab with layers 
perpendicular to the waveguide axis can vary over a wide range, from nearly zero to one, if 
the frequency u of the irradiation is close to the Josephson plasma frequency uj. Meanwhile, 
in the case of layers parallel to the waveguide axis, the transmittance can change significantly 
at frequencies far from uj. 

The paper is organized as follows. In the next section, we discuss the configuration where 
waves propagate in a waveguide along the superconducting layers. We start with the geome- 
try of the problem and present the main equations for the electromagnetic fields both in the 
vacuum regions of the waveguide and also in the slab of the layered superconductor. The 
electrodynamics of layered superconductors is described by nonlinear coupled sine-Gordon 
equation d 3 * 17 ^ . The nonlinearity originates from the nonlinear relation J oc sin ip between 
the Josephson interlayer current J and the gauge-invariant interlayer phase difference ip of 
the order parameter. We emphasize that the nonlinearity can play a crucial role in the prop- 
agation of JPWs even for small wave amplitudes, \ip\ <C 1, when simp can be expanded into 
a series: simp rs ip — p 3 /6. Then we express the transmittance T in terms of the amplitude 
of the incident wave and analyze this dependence for different frequency ranges. In the third 
section, we discuss the same points for the second configuration, when waves propagate in 
a waveguide across the superconducting layers. 

II. PROPAGATION OF WAVES ALONG THE SUPERCONDUCTING LAYERS 
A. Geometry of the problem 

Consider a waveguide of lateral sizes Li and L 2 with a slab of layered superconductor 
of length D inside it. The coordinate system is chosen in such a way that the crystallo- 
graphic ab-plane of the layered superconductor coincides with the xy-plane, and the c-axis 
is along the z-axis. We start our consideration from the case when the electromagnetic 
wave of frequency u propagates in the waveguide along the x-axis, which is parallel to the 
superconducting layers (see Fig. [IJ. 

We now consider an extraordinary incident wave with the magnetic field parallel to the 
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superconducting layers, 

E = {0,0,E z }, H — {H x , H y , 0}. 



(1) 



For simplicity, we assume here that the electromagnetic field is uniform in the ^-direction 
(across the layers) and the electric field is perpendicular to the superconducting layers. Oth- 
erwise, the transformation of the polarization of the wave would occur after its reflection and 
transmission through the anisotropic superconductor. The incident wave is partly reflected 
and partly transmitted through the slab. 




FIG. 1: (Color online) Schematic geometry for JPWs propagating in a waveguide along the su- 
perconducting layers. Note that here S and I stand for superconducting and insulator layers, 
respectively. The green translucent layer (cut-off to show the sample inside) represents the walls of 
the waveguide. 

B. Electromagnetic field in the vacuum regions 

The waveguide has two vacuum regions (see Fig. [I]). In the first one (at x < 0), the 
electromagnetic field can be represented as a sum of incident and reflected waves with 
amplitudes Hi and H r , respectively. Using the Maxwell equations and boundary conditions 
(the equality of the tangential components of the electric field to zero on the waveguide 



5 



walls), one can derive the electric and magnetic field components, 

E { z vl) = [Hi cos(k x x - ujt) 
—H r cos{k x x + tot + a)] sin(k y y), 



j-Hr sm(k x x + cut + a)} cos(k y y), (2) 



k'x 
~k 

+H r cos(k x x + out + a)] sm(k y y). 



= -^[H lC os(k x x-ujt) 



Here 



ky — j , k x — {k ky} , k — , (3) 
lj\ c 

n is a positive integer number that defines the propagating modes in the waveguide, a is the 
phase shift of the reflected wave, and c is the speed of light. 

There is only a transmitted wave of amplitude H t in the second vacuum region (at x > D). 
Here the electromagnetic field components are 

E[f ] = H t cos[k x (x - D) - out + 0\ sm(k y y), 
H (v2) = _k^ Rt sin [ kx ( x _ D )_ ujt + ^ cos (k y y), (4) 
H (v2) = * Ht cos [ kx ( x _ D )_ ut + 0\ sin(^y) 

K 

where (5 is the phase shift of the transmitted wave. 

C. Electromagnetic field in the slab of layered superconductor 

The electromagnetic field in a slab of layered superconductor is related to the distribution 
of the gauge-invariant interlayer phase difference <p(x, y, z, t) of the order parameter (see, e.g., 
Ref. [3]). For the extraordinary wave, this relation is defined by the equations, 
F s K b d 2 H y X 2 ab d 2 H x I dip 



c dtdz ' v c dtdz ' ' : ojjy/e dt 
dH' y dH° x _ U 



simp + 



1 d 2 p 



(5) 



dx dy X 

Here "Ho = 3>o/27r<iA c , <3>o = nch/e is the magnetic flux quantum, \ a b and A c = c/ujjE 1 ! 2 are 
the London penetration depths across and along the layers, respectively, ouj = {Sued J c / he) 1 / 2 
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is the Josephson plasma frequency, J c is the maximal value of the Josephson current density, 
and e is the elementary charge. We do not take into account the relaxation terms because 
they are small at low temperatures and do not play an essential role in the phenomena 
considered here. 

The phase difference is defined by a solution of a set of coupled sine-Gordon equations 
that, in the continuous limit, has the following form (see, e.g., Ref. |3] and references therein): 

x2 d 2 \ r i <9V 1 l2 (&v 9 2 <f\ 

Kbw-^) -2^+sm<? ~K + =0. (6) 



dz 2 



to, at 2 



dx 2 dy 2 



Note that the ^-component of the electric field causes the breakdown of electro-neutrality 
of superconducting layers. This results in an additional, so-called capacitive, interlayer 
coupling. This coupling can significantly affect the properties of the longitudinal JPWs with 
wave vectors oriented across the layers. The dispersion equation for linear plane JPWs with 
account of the capacitive coupling was obtained in Ref. 22, According to this dispersion 
equation, the capacitive coupling is important if the component k z (normal to the layers) is 
close to k — lo/c. For the cases when [k 2 + k 2 ) 1 ^ 2 = co/c, k z = (considered in this section) 
or (k 2 + k 2 ) 1 / 2 ~ k z ~ oj / 'c (considered in the next section), the capacitive coupling can be 
safely neglected due to the smallness of the parameter a = R 2 D e/ sd. Here Rd is the Debye 
length for a charge in a superconductor. 

For the mode in Eq. Q with an electric field E s z proportional to sm.{k y y) and uniform in 
the ^-direction, Eq. ^ is rewritten as 

Equation ([7]) shows that a linear Josephson plasma wave can propagate along the waveguide 
if its frequency exceeds the cutoff frequency u cut , 

u cut = cuj(l + k 2 y \ 2 c y/ 2 . (8) 

Here we consider the weakly nonlinear waves when the Josephson current density J c simp 
can be replaced approximately by J c ((p — y? 3 /6)- We assume that the wave frequency u is 
close to u cu t. In this case, the linear terms in brackets in Eq. ([7]) nearly cancel each other. 
Therefore, in spite of the weakness of the nonlinearity, the term (p 3 plays a very important 
role in the wave propagation. Moreover, one can neglect the generation of higher harmonics 
when the frequency oj is close to u cut (similarly to the effects considered in Refs. [T3fT5]) . 



We seek a solution of Eq. ^ in the form of a wave with the x-dependent amplitude a 
and phase 77, 

ip(x,y,t) = a(x)/j,Ksm[r](x) — cot] sin (k y y), (9) 

where 

v = \y/2,K= |fi 2 - fiLl 1/2 , n = — , n cut = ^. (io) 

Introducing the dimensionless coordinate and the normalized length of the sample, 

x D 

£=—K, = —K, (11) 

A c A c 

and substituting the phase difference <p given by Eq. (|9]) into Eq. ([7]), we obtain two differ- 
ential equations for the functions r)(£) and a(£), 

'/(?) = ~ (12) 

= - a 3 (f) + (13) 

Here a = sign(fi — fi cu t), L is an integration constant, and the prime denotes derivation over 

e 

Using Eqs. ^ and (|9]), we express the tangential components of the electric and magnetic 
fields via the functions a(£) and rj(£), 

E z = —TioXcfinka^) cos[?7(£) — out] sm(k y y), 

H y = Hqiik 2 [a(£) sin[?y(£) - ut}}' sin(k y y). (14) 



Now we can find the amplitudes of the reflected and transmitted waves by solving Eqs. ( 12) 



(13) and matching the tangential components of the electric and magnetic fields at both 



interfaces (at x — and x = D) between the vacuum regions and the layered-superconductor. 
D. Transmittance of the superconducting slab 

In this subsection, we analyze the transmittance T of the finite-length slab of layered 



superconductor in the waveguide. Matching the fields given by Eqs. (14) for the supercon 



ductor and Eqs. ([2]), (pi) for the vacuum regions, we obtain the following three equations for 
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the amplitudes a(0), a(S) and their derivatives on both edges of the layered superconductor 
as well as the expression for the amplitude of the transmitted wave: 

" TL 



[IV 



+ 



+ a(0) 

-TL, 
= 0, 



Here 



h r . 



_a(0) 
a 2 (5) ■ 

a' (5) 
h 2 t =a\5) 



Ahj. 



m 



i, t 



(15) 

(16) 
(17) 
(18) 

(19) 



T-LokXc/jK' 

are the normalized amplitudes of the incident (m = i) and transmitted (m = t) waves, and 
T = K/(k x \ c ). 



Equations (15-17) together with Eqs. (12) and (13) determine the integration constant 



L as a function of the amplitude hi of the incident wave. Using Eqs. (16) and (18), we find 



that the constant L defines directly the transmittance T of the superconducting slab: 

T h? t a 2 (5) T T 



(20) 



hi h 2 hi ■ 

As we show below, the dependence of the transmittance T on the amplitude of the incident 



wave is multivalued due to the nonlinearity of Eqs. (12) and (13). Below, we analyze this 
dependence for the cases of negative and positive frequency detunings (Q — f2 cu t) (negative 
and positive a). 



1. Transmittance of a superconducting slab for negative frequency detunings (a = —1) 
We start our consideration from the case of negative frequency detunings when 

n < n cut . (21) 

As was mentioned above, the linear Josephson plasma waves cannot propagate along the 
superconducting layers under such conditions, i.e., the transmittance of the slab of layered 
superconductor is exponentially small because of the skin effect. However, the nonlinearity 
promotes the wave propagation due to the effective decrease of the cutoff frequency. 



Using Eqs. (12), (13) and conditions rtl5l4l7|, we find the integration constant L and the 



transmittance T given by Eq. (20). Figure 2 shows the numerically-calculated dependence 



T(hi). We consider also the spatial distribution of the amplitude a(£) and the phase tra- 
jectories a' (a). The curves a' (a) are shown in the bottom panels of Figs. M and 111 The 



movement along the spatial coordinate £ (proportional to x, as defined in Eq. (11)) from 
zero to 5 is shown by arrows. The upper panels in Figs. [3] and [4] demonstrate the 3D curves 
for a and a' dependences on the coordinate £. Each trajectory corresponds to some value of 



the amplitude hi. According to Eqs. (16) and (20), the value a(5) defines the transmittance 
of the slab. 

Consider in detail the dependence of the transmittance T on the amplitude hi of the 
incident wave shown in Fig. [2} When increasing the amplitude hi, the transmittance grows 
following the red lower arrows and passing the points 1 and 2 in the T{hi) plot. For such 
amplitudes, the phase trajectories a' (a) represent non-closed curves (see curves 1 and 2 in 
Fig. [3]). At hi « 1.7 the phase trajectory turns into a closed loop (curve 3 in Fig. [3j, 
and transmittance achieves its maximum value T = 1. Then, when increasing hi, the 
transmittance oscillates, achieving the maximum value T = 1 in points 4 and 7 (in Fig. [2j . 
The interesting thing is that there exists a special value of hi ~ 2.4 where the phase trajectory 
a' (a) shrinks into a point (see point 4 in Fig. [3]). This means that, for such value of hi, the 



amplitude a of the electromagnetic wave, given by Eq. (14), is constant inside the slab 
and the phase 77 changes linearly with £, similarly to linear waves. The slab is completely 
transparent in this case (see point 4 in Fig. [2]). 

The continuous change of the transmittance with an increase of hi is abruptly terminated 
in the point 9 (at hi « 7.7) in Fig. [2] Further increase of the amplitude hi results in a jump 
to the upper branch of the T(hi) dependence (to the point 11). Similar jumps occur again 
and again when further increasing the amplitude hi of the incident wave. 

Let us now consider the T{hi) dependence when decreasing the incident wave amplitude 
hi, following the blue upper arrows in Fig. [2] The decrease of hi causes an increase of the 
transmittance T untill we achieve the situation when the slab is completely transparent. Here 
a branch of the T(hi) dependence breaks and further movement along the curve is possible 
only after a jump to the lower branch. As one can see from Fig. [2j an additional decrease of 
the amplitude hi results in a similar behavior until we reach the point 10 (hi ~ 6.8), where 
the last jump to the lower branch (to the point 8) occurs. Further decrease does not exhibit 
additional jumps. 

Thus, for frequencies Q less than the cutoff frequency fi cn t, the transmittance of a super- 
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FIG. 2: (Color online) The transmittance T versus the normalized amplitude hi of the incident 
wave for negative frequency detuning, £7 — fj cut < 0. Movement along the lower red (upper blue) 
arrows shows the evolution of the transmittance when increasing (decreasing) the amplitude hi. 
The numbers near the points on the T{hi) curve correspond to the same numbers of the phase 
trajectories a' (a) shown in Figs. [3] and [4] The inset demonstrates the enlarged region near the first 
three maximums in the T(h{) dependence. The values of the parameters are: k = 0.1, 5 = 1.5 
(D = 15A C ), A c = 4 • 10" 3 cm, X ab = 2000 A, lo^tt = 0.3 THz, m = 1, and L x = 0.1 cm. 
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FIG. 3: (Color online) The phase trajectories a! {a) for negative frequency detuning, Q — f2 cut < 0. 
These trajectories correspond to points from 1 to 6 in the T{hj) plot shown in Fig. [2] The upper 
panel demonstrates in 3D the dependences of a and a' on the spatial coordinate £ along the slab 
(for points 1, 3, and 6 in Fig. ^1). The parameters here are the same as in Fig. ^1 



conducting slab has a hysteretic dependence with jumps on the amplitude of the incident 
wave. 

In Fig. [5j we illustrate the spatial distribution of the electric field amplitude inside the 
waveguide for three values of the incident amplitude h^. Panel [51(a) corresponds to the point 
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FIG. 4: (Color online) The same as in Fig. [3] for the points 7-12 shown also in Fig. [2] 



10 in Fig. |2j where the layered superconductor is almost transparent. This means that the 
reflected wave is almost absent and the interference pattern in the vacuum regions is too 
blurry. The panels [5](b) and[5|c) correspond to the points 11 and 12 in Fig. [2j They show a 
pronounced interference pattern for the region x < 0. Note that the field distribution along 
the ?/-axis corresponds to n\ = 1. 
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FIG. 5: (Color online) The spatial distribution of the electric field amplitude (normalized to the 
incident wave amplitude Hi) inside the waveguide. Panels (a), (b), and (c) correspond to the points 
10, 11, and 12 in Fig. |2j respectively. The color determines the value of the amplitude. The vertical 
straight lines show the edges of the superconducting sample. The parameters are the same as in 
Fig. [2} 
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2. Transmittance of a superconducting slab for positive frequency detunings (a = +1) 



Now we discuss the case of positive frequency detunings when 

n > n cut . (22) 

Under this condition, contrary to the case Q < fi cut , even the linear Josephon plasma 
waves can propagate in the superconductor along the layers. Therefore, the transmittance 
is not exponentially small in the linear regime. It vary over a wide range, from zero to one, 
depending on the relation between the length of the slab and the wavelength. Figure [6] shows 
the dependence of the transmittance on the incident wave amplitude for different frequency 
detunings. For proper comparison of the curves, we use here another normalization of the 



incident wave amplitude, independent of k and, therefore, different from Eq. (19). Namely, 
we plot all the curves as functions of the parameter hi ■ k. One can see that the transmittance 
T depends strongly on k and a for small values of hi ■ K. However, all the dependences of T on 
the incident wave amplitude show similar behavior at larger hi ■ K. This means that, for high 
enough amplitudes Hi, the nonlinearity plays the most important role in the transmissivity, 
regardless of the frequency detuning. 



E. Mechanical analogy 



The problem discussed here has an interesting mechanical analogy. Indeed, Eqs. (12) 



and (13) coincide with equations describing a motion of a particle of mass one in a centrally 



symmetric field. The amplitude a(£), the phase and the coordinate £ along the layers of 
the superconductor play the roles of the radial coordinate, polar angle, and time, respectively. 



The constant L in Eqs. (12) and (13) can be considered as the conserved angular momentum 
of the particle. 



Integrating Eq. (13), we find the energy conservation law, 

^ + U cS {a) = £, 



(23) 



with the effective potential energy 

U eS (a) 



L 2 a 2 a 4 

2c7 2+a Y + T 



(24) 
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FIG. 6: (Color online) Dependence of the transmittance T on the value of hi ■ K for different fre- 



quency detunings. According to Eq. (19), the argument hi -K is proportional to the non-normalized 



amplitude Hi of the incident wave and is independent of the frequency detuning. The red, or- 
ange, and yellow curves (the upper curves at small amplitudes) correspond to positive frequency 
detunings with k = 0.2, 0.1, and 0.05, respectively. The purple, blue, and green curves (the lower 
curves at small amplitudes) correspond to negative frequency detunings with k = 0.2, 0.1, and 0.05, 
respectively. Other parameters are the same as in Fig. [2] 



The first term in Eq. (23) is the kinetic energy of the radial motion, £ is the total energy, the 



first term in Eq. (24) is the centrifugal energy, and the last two terms in Eq. (24) describe 



the potential of the central field. 
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FIG. 7: (Color online) Dependence of the potential C/ e ff given by Eq. (24) on the radial coordinate a. 



The movement of the particle in this potential is the mechanical analog for the spatial distribution 
of the amplitude a of the electric field in the superconductor. The solid red curve corresponds to 
the positive frequency detuning and the dashed blue curve corresponds to the negative one. The 
value of the constant L is 0.1, and the other parameters are the same as in Fig. [2] 



The dependence U e g(a) illustrated in Fig. [7] shows that the potential energy is a single- 
valued function of a, although the dependence T{hi) is multivalued (see Fig. [2]). This feature 
seems to be paradoxical, because the particle motion in any potential field is unambiguously 
denned by the initial conditions. However, an assignment of the value of hi of the incident 



wave amplitude in Eqs. (15 -17) not always mean an imposition of definite initial conditions 
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for the particle motion. To demonstrate this, let us consider the inverse problem. We seek 
the amplitude hi of the incident wave that is necessary to obtain a definite amplitude h t of the 



transmitted wave. From Eqs. ( 16 ) and ( 18 ) , we see that the value of h t defines unambiguously 



the angular momentum L = hq/T of the particle. According to the motion equations (12) 



(13) and the boundary conditions Eqs. (15-17), the dependence of hi on the amplitude ht 



of the transmitted wave and, correspondingly, the dependence of the transmittance 



T = h\lh\ 



on h t are single- valued (see Fig. |8J). Nevertheless, the dependence h{(ht) is nonmonotonic. 
Therefore, the dependence T(hi) appears to be multiple-valued. 



III. PROPAGATION OF WAVES ACROSS SUPERCONDUCTING LAYERS 

Consider now the case when the electromagnetic wave propagates in the waveguide across 
the superconducting layers, i.e., the crystallographic c-axis of the superconductor and the 
z-axis are parallel to the waveguide axis (see Fig. Km. Let the incident extraordinary wave 
have the following polarization: 

E = {E x ,E y ,E z }, H = {H x ,H y ,0}. (25) 



In the first vacuum region (at z < 0) the electromagnetic field can be represented as a 
sum of incident and reflected waves with amplitudes Hi and H r , respectively. There is only 
the transmitted wave of amplitude H t in the second vacuum region (at z > D). Using the 
Maxwell equations and boundary conditions (the equality of the tangential components of 
the electric field to zero on the waveguide walls), one can derive the electric and magnetic 
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Transmitted amplitude h t 

FIG. 8: (Color online) Solution of the inverse problem (the dependence of the incident wave am- 
plitude hi and the transmittance T = h^/hf on the transmitted wave amplitude ht). The values of 
the parameters and the numbers of the points on the curves correspond to the numbers in Fig. [2] 
The curves show the nonmonotonic behavior of hi(ht)- This results in the multivalued dependence 
T(hi). 
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FIG. 9: (Color online) Schematic geometry for JPWs propagating in a waveguide across the su- 
perconducting layers. Note that here S and I stand for superconducting and insulator layers, 
respectively. The green translucent layer (cut-off to show the sample inside) represents the walls of 
the waveguide. 



field components in the first and second vacuum regions: 

[Hi sin(k z z — ut) — H r sin(&; 2 ,2 + ut + a)} 



E (vi) _ k xk 



k 2 



x cos(k x x) sm(k y y), 



k k 

Ej. vl > = [Hi sin(k z z - ut) - H r sin(k z z + ut + a)} 
k 2 



x sin(k x x) cos(k y y), 



k,. 



H^ l) = —jL[Hi sia(k z z - ut) + H r sin(k z z + ut + a)} 



x sm(k x x) cos(k y y), 



(vi) 



H y vV} = [H sm(k z z - ut) + H r sm(k z z + ut + a)} 



(26) 



E (v2) 



\\i >r< /u 9 



H t sm.[k z (z — D) — ut + (3] cos(k x x) sxn(k y y), 



E y v2) = -j^H t sm[k z (z -D)-ut + f3] sm{k x x) cos{k y y), 



Ui.v2) 



k, 



H t sin[k z (z — D) — ut + (3] sm(k x x) cos(k y y), 



A r 



H (v2) = x Rt sin [ kz ( z _ D )_ ujt + p] cos (k xX ) sm(k y y). 



(27) 
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Here 

k x = — , k y = — , k z = [k 2 — k 2 — kTj , (28) 

rii and n 2 are nonnegative integer numbers that define propagating modes in the waveguide 
(they cannot be equal to zero simultaneously), a and /3 are the phase shifts of the reflected 
and transmitted waves. We do not write the equation for the z component E z of the electric 
field since it is not used further in the text. 

For the region occupied by the layered superconductor, we seek a solution of Eq. ^ in 
the form of a mode, 

<p(x,y,z,t) = a(z)\l - Vt 2 \ l/2 

x sm(k x x) sin(k y y) sin[r)(z) — cot], (29) 

with the z-dependent amplitude a(z) and phase r)(z). 

We introduce the dimensionless coordinate ( and the normalized length of the sample as 

C=t~i 5=—, (30) 

Aab *ab 

where 

- = A C ^W h = W + K) V2 ■ (31) 
Substituting the phase difference ip in the form of Eq. (29) into Eq. ([6]), we obtain two 
differential equations for the phase r](() and the amplitude a((), 

'/(C) = -J^y (32) 

=<•«) + ^ + (33) 
where L is an integration constant, prime denotes derivation over £, and 

h(C) = -sign(fi - l)a(C) - ^« 3 (0- (34) 

Note that the cutoff frequency for linear waves propagating across the superconducting layers 
coincides with the Josephson plasma frequency ujj. Therefore, for the geometry considered 
in this section, the transmittance of the superconducting slab is sensitive to the frequency 
detuning (fl — 1). 
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Equations Q, g, and Q) allow one to express the electromagnetic field inside the slab 
via the functions h(Q and 77(C)- The tangential components of the electric and magnetic 
fields are: 

ei = -n v 1 

k k\\K 

x {h(Q cos[?7(C) — urt]} cos(k x x) sxn(k y y), 

E s = Y ^ y k z |1 - £l 2 \ 
y k k\\R 

x {h(() cos[?7(C) — tut]} sin(k x x) cos(k y y), (35) 
II - ft 21 



H^, — Hoky- 



k\\K 



MO 



x sin [77(C) — ojt] sin^a;) cos(k y y), 

y K\\K 

x sin[r/(C) — out] cos(k x x) sxn(k y y) 



with 



k 2 \ 



ab - 



k? 



-K. 



(36) 



Matching the fields given by Eqs. (35) for the superconductor and Eqs. (26), (27) for the 



vacuum regions, we obtain the following three equations for the amplitudes a(0), a(5) and 
their derivatives on both edges of the layered superconductor, as well as the expression for 
the amplitude of the transmitted wave: 

2 



r 2 [h 



+ 



TL 



h(0) + MO) 

h 2 (5) = TL, 
h'{5) = 0, 
h 2 t = h 2 (5). 



(37) 

(38) 
(39) 
(40) 



Here 



H„ 



k\\K 



H n k\l~m 



m = i,t 



(41) 



are the normalized amplitudes of the incident (m = i) and transmitted (m = t) waves. 



These equations, together with Eqs. (32) and (33), determine the integration constant L as 
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a function of the amplitude hi of the incident wave. Similarly to the previous section, the 
constant L defines directly the transmittance T of the superconducting slab, 



Hi h 



L. (42) 



Figure 10 shows the numerically-calculated hysteretic dependence of the transmittance 



T on the amplitude hi of the incident wave for the case of negative detuning, Q — 1 < 0. 



The red curve close to the abscissa in Fig. 10 shows the low- amplitude branch of the T(hi) 



dependence. The low-amplitude waves cannot propagate across the superconducting layers 
at Q < 1 and, therefore, the transmittance is exponentially small on this branch. When 
increasing the amplitude hi of the incident wave, the jump from the low-amplitude branch 
of the T(hj) dependence to the high-amplitude branch occurs at hi = (128/243) 1 / 2 . The high- 



amplitude branch is shown by the blue solid curve in Fig. 10 The high- amplitude solutions 
describe nonlinear Josephson plasma waves that can propagate in the layered superconductor 
for even negative frequency detunings. Changing the amplitude hi one can control the 
relation between the wavelength and the length of the slab, and the transmittance varies 
from nearly zero to one depending on this relation. Thus, one can obtain total transparency 
of the slab choosing the optimal value of the amplitude hi. 



Figure 11 shows the T{hi) dependence for the case of positive detunings, Vt — 1 > 0. In this 
case, even linear JPWs can propagate in the superconductor across the layers. Therefore, 
the linear transmittance can vary from nearly zero to one depending on the relation between 



the length of the slab and the wavelength. The analysis of Eqs. (32-34), (37-40), and (42) 
shows that the dependence T{hj) is single-valued if the frequency detuning exceeds some 
threshold value. For small T, this threshold value is defined by the following asymptotic 
equation: 

/ D\ c kn \ 2 

^-^b^) • (43) 



Such a reversible T{hi) dependence is shown by the dotted curve in Fig. 11 The hysteresis 
in the T(hj) dependence exists for frequencies smaller than the threshold value, for Q < Othr- 
In this case, the complete transparency of the sample can be observed if the incident wave 
amplitude hi is first increased and then decreased. 



Note that plots in Figs. 10 and 11 are very similar to T{hi) curves obtained in Ref. HE] 
for the case of an infinite slab. It is not surprising because the statements of these two 
problems differ in not very important details. Namely, the nonlinear waves running in the 
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Incident amplitude hi 



FIG. 10: (Color online) Dependence of the transmittance T on the amplitude h% of the incident wave 
for negative frequency detuning, Q — 1 = —5 • 10~ 5 . Arrows show the change of the transmittance 
when changing hi. The inset shows the enlarged region near the points 1 and 2. The values of the 
parameters used here are: 5 = 2, A c = 4 • 10 -3 cm, A a & = 2000 A, ujj/2tt = 0.3 THz, n\ = n-i = 1, 
and L\ = L<i = 0.1 cm. 
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0.01 0.02 
Incident amplitude /?; 



FIG. 11: (Color online) Dependence of the transmittance T on the amplitude hi of the incident wave 
for different positive values of the frequency detuning: Q — 1 = 4.97 • 1CP 3 , or 6 /it = 1.21 (dotted 
curve); Q - 1 = 4.74 • 1CT 3 , or 5/tt = 1.24 (solid curve); O - 1 = 1.64 • 1CT 3 , or Ij-K = 2.11 (inset). 
Arrows show the change of the transmittance when changing the amplitude hi of the incident wave. 
The sample length is D = 4.3 • 1CP 5 cm. Other parameters are the same as in Fig. 10 
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infinite slab at some angle with respect to the superconducting layers were considered in 
Ref. Uni whereas the nonlinear waves standing in the x- and y-directions and propagating 
only across the layers are studied in this section. Therefore, the sets of equations for finding 
the transmittance in these two problems differ in the numerical coefficients only. 

IV. CONCLUSION 

In this paper, we have studied theoretically the phenomenon of self-induced transparency 
of finite-length slabs of layered superconductors placed in a waveguide with ideal walls. 
This phenomenon can be observed for two geometries: when the superconducting layers are 
parallel to the waveguide axis or perpendicular to them. We show that the transmittance 
of a superconducting slab is very sensitive to the amplitude of the incident wave due to the 
nonlinear dependence of the Josephson c-axis current on the interlayer phase difference of the 
order parameter. A very interesting feature of the predicted phenomenon is the hysteretic 
behavior of the T(hj) dependence. It is important to note that the tunable transmittance can 
vary over a wide range, from nearly zero to one, when changing the incident wave amplitude, 
even in the case of weak nonlinearity when the interlayer phase difference is small. 
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